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Abstract 

Oscillations in the transcriptional activator NF-kB localized in the nucleus have been 
observed when a cell is stimulated by an external agent. A negative feedback based on 
the protein IkB whose expression is controlled by NF-kB is known to be responsible 
for these oscillations. We study NF-kB oscillations, which have been observed both for 
cell populations by Hoffmann et al. (2002) and for single cells by Nelson et al. (2004). In 
order to study cell-to-cell variability we use Gillespie's algorithm, applied to a simplified 
version of the model proposed by Hoffmann et al. (2002). We consider the amounts of 
cellular NF-kB and activated IKK as external parameters. When these are fixed, we 
show that intrinsic fluctuations are small in a model with strong transcription, as is the 
case of the Hoffmann et al. (2002) model, whether transcription is quadratic or linear in 
the number of NF-kB molecules. Intrinsic fluctuations can however be large when tran- 
scription is weak, as we illustrate in a model variant. The effect of extrinsic fluctuations 
can be significant: cell-to-cell fluctuations of the initial amount of cellular NF-kB affect 
mainly the amplitude of nuclear NF-kB oscillations, at least when transcription is linear 
in the number of NF-kB molecules, while fluctuations in the amount of activated IKK 
affect both their amplitude and period, whatever the mode of transcription. In this case 
model results are in qualitative agreement with the considerable cell-to-cell variability of 
NF-kB oscillations observed by Nelson et al. (2004). 
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1 Introduction 



For a number of biological processes cell-to-cell variability can be very large, even for clonal 
populations, as pointed out by Dower and Qwarnstrom (Dower and Qwarnstrom, 2003), who 
studied immune response signaling systems connected to the family of TIR (Toll/interleukin- 
1) receptors. In such cases average cell behavior is a poor indicator of the way individual cells 
respond. When for example cell expression as a function of inducer or enhancer is stochastic, 
graded average response can hide all or nothing single cell response, as shown in a number 
of studies (Fiering et al., 1990, 2000; Hume, 2000; Ko, 1992; Nair et al., 2004). It is there- 
fore important that in both experiment and modeling cell-to-cell variability be investigated. 
In particular, one needs to be attentive to this fact when modeling a biological system with 
a set of chemical rate equations for concentrations: these by definition encapsulate average 
behavior only. A case in point is provided by recent experiments and modeling of NF-kB 
activation as a result of TNFa stimulation of human Jurkat T cells, U937 monocytes, and 
mouse fibroblasts (Hoffmann et al., 2002), or human cervical carcinoma and type-S neurob- 
lastoma cells (Nelson et al., 2004). Hoffmann et al. (2002) measure average cell behavior 
and model it with a set of chemical rate equations for concentrations, whereas Nelson et al. 
(2004) examine the same model but give experimental results on single cells only, some of 
which show large cell-to-cell variations. It is difficult to reconcile well defined average oscil- 
latory behavior with cell-to-cell variability in both amplitude and period of oscillations. Our 
aim is to understand (within the confines of the model of Hoffmann et al. (2002)) the origin 
of the observed (Nelson et al., 2004) cell-to-cell fluctuations of nuclear NF-kB oscillations. 

The fluctuations which underlie cell-to-cell variability can be separated into intrinsic and 
extrinsic (Elowitz et al., 2002; Raser and O'Shea, 2004). Intrinsic fluctuations are those 
which result from random occurrence of a defined set of reactions at random times, leading 
to cell-to-cell variability for genetically identical cells in identical, fixed environments. Ex- 
trinsic fluctuations can have multiple origins, such as fluctuations in cytoplasmic and nuclear 
volumes, in the initial state of the cell, signaling cascade fluctuations, or variations in the 
number of regulatory molecules. 

We should point out that NF-kB oscillations are a transient phenomenon, unlike the 
oscillatory process involved in circadian rhythms. Under persistent stimulation of the cell, 
activation of IKK itself is temporary, either through adaptation (Hoffmann et al., 2002) or 
through the action of inhibitory factor A20 (Lipniacki et al., 2004). Thus mathematically the 
cell which is in a steady state before stimulation returns to a steady state some time (many 
hours) after stimulation. The amplitude and period of oscillations change as time progresses, 
and they disappear ultimately. 

We use stochastic simulations (Gillespie, 1977) of the rate equations and show the follow- 
ing: 
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- intrinsic fluctuations are very small for the model parameters used in the literature (Hoff- 
mann et al., 2002; Nelson et al., 2004), which correspond to strong transcription, and they 
affect only late nuclear NF-kB oscillations. Intrinsic fluctuations can however be large if 
promoter binding is weak, with rate constants for transcription and translation significantly 
different from those used by Hoffmann et al. (2002) 

- extrinsic fluctuations can be significant: fluctuations in the initial amount of activated IKK 
lead to both amplitude and periodicity fluctuations of NF-kB oscillations, which are features 
shown by the single cell data of Nelson et al. (2004). This result, for the model parameter 
set of Hoffmann et al. (2202), is independent of whether transcription is linear or quadratic 
in the amount of nuclear NF-kB. However for fluctuations in the initial amount of cellular 
NF-kB cell-to-cell variability depends on the mode of transcription: only for quadratic tran- 
scription are both amplitude and period affected. For linear transcription only the amplitude 
of oscillations varies. 

These results on cell-to-cell variability make clear the importance of determining exper- 
imentally the form of the transcriptional relationship between NF-kB and IkBoi, and the 
actual values of transcriptional and translational rate constants. 

In the next section we describe the model. This is followed by a section of results. Here 
we discuss intrinsic noise and study a variant of the Hoffmann et al. (2002) model with a 
different transcription mechanism in order to gain more insight into the generation of intrinsic 
noise. We also describe the results of fluctuations in the initial amount of NF-kB, and in the 
number of activated IKK molecules. We conclude with a discussion in which we speculate 
on the functional role of a threshold in the number of NF-kB oscillations, as cell stimulus, 
and therefore the amount of activated IKK decreases in intensity. 

2 Description of model 

NF-kB (nuclear factor kB) plays an important regulatory role in cell immune response and 
survival (for a review, see Tian and Brasier, 2003). Under normal conditions, it is sequestered 
in the cell cytoplasm through binding to IkB. When the cell is stimulated through for exam- 
ple TNFa, the ensuing signaling cascade leads to activation of IKK, which in turn through 
binding to NF-kB and IkB, leads to ubiquitination of IkB, and the freeing of NF-kB. 
NF-kB enters the cell nucleus (Carlotti et al., 1999), and regulates the expression of a num- 
ber of genes, in particular the gene expressing IkB. IkB protein enters the nucleus, binds to 
NF-kB, and that complex is exported. Thus NF-kB which initially was sequestered, then 
released after IKK activation, ends up by being sequestered again. As a result, the concen- 
tration of nuclear NF-kB shoots up shortly after stimulation, and then decreases because of 
the negative feedback associated with IkB expression through NF-kB itself binding to the 
relevant promoter. The process can repeat itself if the time scales associated with nuclear 
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import and export, transcription and translation are right. Indeed measurements show that 
the concentration of nuclear NF-kB oscillates over several hours after the onset of a lasting 
stimulation (Hoffmann et al., 2002; Nelson et al., 2004). 

For the points we wish to make concerning cell-to-cell variability, a pared down version 
of the model of Hoffmann et al. (2002), which only includes the IkB isoform responsible for 
NF-kB oscillations, suffices. The original model includes 64 ordinary differential equations 
for concentrations, with reaction rate constants determined or constrained from experiment, 
and others constrained by fits to the data obtained with the use of GEPASI software ( Mendes, 
1993). 

There are three isoforms of IkB, the main one being IkBol, whose transcription depends 
strongly on the presence of nuclear NF-kB. For the other two, transcription is not regulated 
by NF-kB; Hoffmann et al. (2002) were interested in elucidating the role of these other two 
isoforms in damping NF-kB oscillations. Therefore their system of equations is extensive, 
each isoform participating in the same set of reactions (except for NF-kB dependent tran- 
scription), but without any interaction between isoforms themselves. In our simplified version 
(the number of reactions is about 1/3 of those of the original model), we will only consider 
the dominant isoform InBa, which is solely responsible for the occurrence of NF-kB oscil- 
lations, and the only one whose knockout is lethal (Gerondakis et al., 1999). 

Since we are interested in intrinsic noise issues, a model based on ordinary differential 
chemical rate equations is inadequate. We use Gillespie's algorithm (Gillespie, 1977, 2001) 
to study the impact of including both the random occurrence of chemical reactions, and fluc- 
tuations in external parameters on the shape of NF-kB oscillations. 

We should mention another model of NF-kB oscillations (Lipniacki et al., 2004), which 
has as its starting point the Hoffmann et al. (2002) model, but focuses in addition on the role of 
A20, which is regulated by NF-kB, and acts as an inhibitor of IKK, itself activated as a result 
of cell stimulation by an external agent. There is also a discussion of NF-kB oscillations by 
Monk (2003) in the context of a model of ordinary differential equations with time delays. 

The reactions we have used in our model are summarized in Table 1 together with the 
relevant rate constants, which correspond to the choices made by Nelson et al. (2004, sup- 
plementary material)) for the reactions of Hoffmann et al. (2002). Transcription is linear in 
the number of NF-kB molecules (Lipniacki et al., 2004; Nelson et al., 2005), rather than 
quadratic as in the original model (Hoffmann et al., 2002). This change, from quadratic to 
linear, diminishes the sensitivity of NF-kB oscillations to fluctuations in the initial amount 
of NF-kB present in a cell (Barken et al., 2005; Nelson et al., 2005). We will give results for 
both choices of transcription in the next section. The model as it stands still has some unsat- 
isfactory features, namely the amount of transcribed molecules rises linearly (or quadratically 
as in Hoffmann et al. (2002)) with the number of NF-kB molecules entering the cell nu- 
cleus, without any mechanism of saturation, and the number of iKBa mRNA molecules is 
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excessive, of the order of a thousand in the steady state (for the simulation underlying Figure 
1), rising to much higher numbers (several times ten thousand) once IKK is activated. We 
have however checked, as we will discuss in section 3.1., that as a result of decreased tran- 
scription and a corresponding increased translation, the number of IkBcx mRNA molecules is 
reduced, while NF-kB oscillations themselves are not affected, because the amount of NF- 
kB present in the cell is large; the intrinsic noise in the system, however, becomes larger, the 
more so as the transcription rate decreases. 

In the model, NF-kB is not degraded, and therefore the number of NF-kB molecules 
remains constant through the formation and disintegration of complexes in which it partici- 
pates. We use cells with volumes of 1000 /i 3 , divided equally between cytoplasm and nucleus. 
We take an amount of 0.05yuM relative to the total cell volume for NF-kB. Thus the number 
of NF-kB molecules initially introduced, corresponding to 0.05/iM for the whole cell, is 
30,000. After an equilibration period of 3000 minutes, we follow Hoffmann et al. (2002) in 
introducing an equal amount of IKK (0.05 fiM), which represents activation of the signaling 
pathway through cell exposure to TNFa. 

The shape of NF-kB oscillations will for a given cell volume depend on (all other pa- 
rameters being fixed) the fraction of volume occupied by the nucleus (Lipniacki et al., 2004), 
because reactions take place in both cytoplasm and nucleus and there is nuclear import and 
export. Their shape depends as well on the comparative amount of initial NF-kB and IKK 
activated after cell stimulation. 

NF-kB oscillations, which are more or less damped, are sensitive to a number of model 
parameters such as the rate of transcription of IkBcx or its translation, or IkBo: mRNA degra- 
dation, which are unknown and determined from fits to experiments. Thus it is not always 
possible to evaluate the significance of the precise form of the oscillations from comparing 
model results with experiment. 

3 Results. 

Data are obtained in different cell lines (Hoffmann et al., 2002; Nelson et al., 2004) with 
different techniques and therefore comparisons can be hazardous. Nevertheless, the data on 
NF-kB signaling show well defined oscillations ( of period close to 100 minutes) when av- 
eraged over cells, when only isoform IkBo is present (Hoffmann et al., 2002, supplementary 
material, fig. SI). However the single cell data of Nelson et al. (2004) show large cell-to-cell 
variations in both amplitude and periodicity, which if typical, would be expected to broaden 
or even wash out average NF-kB oscillations. It is in this context that we will next discuss 
our results. We will examine three possible sources of cell-to-cell variability: intrinsic fluc- 
tuations on one hand, and on the other extrinsic fluctuations in either the initial number of 
NF-kB molecules or of activated IKK molecules. 
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3.1 Intrinsic noise. 



Intrinsic noise is small in the models of Hoffmann et al. (2002) and Nelson et al. (2004) given 
the rate constants they employ. Whether transcription is linear or quadratic in the number 
of NF-kB molecules, intrinsic fluctuations become significant only about 600 minutes after 
IKK activation. In Figure 1 we show, for the case of linear transcription, and for the param- 
eters of Table 1, nuclear NF-kB oscillations as a function of time after activation by IKK 
after 3000 minutes of equilibration. For this figure the initial numbers of NF-kB and acti- 
vated IKK molecules are fixed. Thus any fluctuation is purely intrinsic, due to the random 
order in which reactions occur in any given cell. The first three or four oscillations have a pe- 
riod of about 100 minutes, thereafter the period lengthens. The figure shows average behavior 
and the results for three randomly chosen single cells. The comparison illustrates the point 
that fluctuations are very small: in each cell oscillations for 8 hours after activation are almost 
identical to the average. (The rapid small amplitude fluctuations over times less or much less 
than a minute are washed out in the representation of Figure 1 .) Thus cell-to-cell variability 
is very small. The main reason is that the number of molecules is large: with values of the 
order of 0.05 fxM for initial numbers of NF-kB and activated IKK molecules, and cells 
of volume 1000 fi 3 , their numbers are in the tens of thousands. Moreover transcription, the 
strength of which increases with the number of NF-kB molecules, is strong. If in the models 
considered, one decreases, while keeping all other rate constants unchanged, the rate of tran- 
scription by merely a factor of 2 (for linear transcription) or 4 (for quadratic transcription) , 
NF-kB oscillations beyond the first one are very weak and disappear. 

Nevertheless the strength of intrinsic fluctuations indeed varies with the values of the rate 
constants. As already mentioned in Section 2, we have decreased, within the model with 
linear transcription, the rate of transcription of IkBq. mRNA by a factor of 10, in order to 
have fewer mRNA, and increased translation by the same factor, thus keeping the number 
of iKBa proteins essentially unchanged. This model still yields good oscillations in nuclear 
NF-kB but with increased cell-to-cell variability, as shown in Figure 2, and as compared to 
the single cell curves in Figure 1. The decreased copy number of IkBcx mRNA thus leads to 
noise in its own oscillations that is reflected in the behavior of NF-kB. The intrinsic noise 
is increased (up to factors of two) over the original model. We will investigate further the 
connection between strength of transcription and intrinsic noise in the following section with 
a more appropriate model which includes explicit promoter binding. 

3.1.1 Model variant with DNA. 

The model of Hoffmann et al. (2002) treats transcription in an oversimplified way without 
including DNA. In order to examine the effect of the stochastic binding and unbinding of 
NF-kB to the promoter site of the DNA on the level of intrinsic noise we have considered 
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the following model: the transcription of the mRNA coding for InBa, denoted by Ictr, occurs 
through the reactions. 

DNA + Nn ^ DNA* (1) 
kb 

DNA* ^% DNA + Nn + Ictr (2) 

Since the sum of the concentrations of the DNA and the complex DNA* is conserved, we 
have, c D + c D * = c where c = 2/300000 piM given the volume of the nucleus (500yU 3 ) and 
assuming a diploid cell. If the first reaction above occurs rapidly and is hence, in equilibrium, 
we can combine the two equations in the spirit of Michaelis-Menten to obtain 

Nn tr ^ f Nn + Ictr 

where the effective transcription rate that depends on the nuclear NFkB concentration Nn is 
given by 

c tr' 2 

tr2 > eff ~ %Tmi • 

Since the maximum value of Nn is 0.1/iM we chose kb/kf = 0.2 /iM; changing the value 
up to 1.0 fxM had negligible effect on the results given below. In order to obtain similar 
transcription rates as for the linear model we chose tr' 2 x c to be equal to the value of tr 2 used 
in the original model given in Table 1 . In order to ensure that the binding and unbinding of Nn 
occur sufficiently rapidly we chose kb = 1200 min^ 1 and the translation rate trl is increased 
by a factor of 2 from the value given in Table 1 . With this high rate of promoter binding, one 
is in a situation, as for the models of linear or quadratic transcription (Hoffmann et al., 2002; 
Nelson et al., 2004), where the relevant promoter is strongly activated and transcriptional 
stochasticity is small (Raser and O'Shea, 2004). This is confirmed by the simulation reported 
in Figure 3, which shows the concentration and standard deviation of nuclear NF-kB as a 
function of time after activation by IKK after 3000 minutes of equilibration. The noise is 
small in the first four oscillations and increases at larger times, in the same way as shown in 
Figure 1 for linear transcription with the parameter set of Table 1 . 
Transcriptional noise. 

Moving away from the regime of strong transcription embodied in the models with linear 
or quadratic transcription (Hoffmann et al., 2002; Nelson et al., 2004), one can now explore 
what happens when promoter binding gets weaker, with a rate that gets closer to the rates 
which determine InBa protein ubiquitination, namely rates rl and r4 (cf. Table 1). Tran- 
scriptional fluctuations are then no longer averaged out over the corresponding time scales 
(Lewis, 2003). One can do so since, contrary to most rate constants which are determined or 
constrained by experiment (Hoffmann et al., supp. material, 2002), those of transcription and 
translation, tr2 and trl respectively (cf. Table 1), result from fits to the oscillation data. How- 
ever, as one decreases the rate of promoter binding, it is necessary to increase correspondingly 
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the rate of translation in order to maintain well defined nuclear NF-kB oscillations. Figure 
4 illustrates the impact of weak transcription on the strength of intrinsic noise, for a binding 
rate of kb = YlmirC 1 , a hundredfold decrease as compared to the case of strong transcription 
(see above), and a corresponding 200 fold increase of the translation rate trl of Table 1. The 
intrinsic noise is much more significant than in Figure 3, and, as shown in Figure 5, it affects 
both the amplitude and periodicity of NF-kB oscillations. Indeed, the intrinsic noise is suf- 
ficiently large that the oscillations obtained from averaging the stochastic simulations is more 
damped than those obtained from the deterministic system. This underscores the potential 
pitfalls of deducing rate constants by fitting with deterministic models. 

One could investigate the effects of further reducing promoter binding. For our purpose 
which is to establish and clarify the possible sources of the observed cell-to-cell variability 
(Nelson et al., 2004) of NF-kB oscillations, it is sufficient to emphasize that intrinsic noise 
can be significant, as important as extrinsic noise, if we do not confine ourselves to the fitted 
values of transcription and translation rates of the set of reactions considered by Hoffmann et 
al. (2002). 

3.2 Extrinsic noise. 

Since the data indicate large cell-to-cell fluctuations (Nelson et al., 2004), these must arise 
from fluctuations in external parameters, at least for the models with strong transcription, for 
which, as shown above, intrinsic noise is small. We consider two possible sources here, the 
number of initial NF-kB present in a given cell, which varies between cells in an analogous 
situation (Dower and Qwarnstrom, 2003), and the amount of IKK, which can vary due to 
signaling cascade fluctuations. The models considered are those with transcription quadratic 
in the number of NF-kB molecules (Hoffmann et al., 2002) and its version with linear tran- 
scription (Nelson et al., 2004) with the rate constants of Table 1. For quadratic transcription, 
the linear transcription reaction in Table 1 is replaced by 2Nn — > 2Nn + Ictr with rate con- 
stant trl = 1.027125 / uM~ 1 mm~ 1 . For Figure 7 moreover rate constants trl and tr 3 of Table 
1 are multiplied by 2. 

3.2.1 NF-kB fluctuations. 

As alluded to in section 2, there is a significant difference here between the models with linear 
and quadratic transcription (Barken et al., 2005; Nelson et al., 2005). For linear transcription, 
moderate fluctuations in the amount of NF-KB(see below) affect the amplitude only of single 
cell oscillations, thus maintaining well-defined average oscillations. For quadratic transcrip- 
tion, both amplitude and period of single cells fluctuate, leading to reduced or washed out 
average oscillations. The following discussion explains this point. 

In Figure 6, for linear transcription, are shown results for NF-kB oscillations for the 
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average over many cells for a Gaussian distribution of initial number of NF-kB (with average 
30000, a =5000) (black curve) and single cell behavior when initial NF-kB is one or two 
a away from the average. IKK remains fixed in all cases at IKK=30000. Average NF-kB 
oscillations in Figure 6 are of course very similar to those in Figure 1, diverging from them 
only beyond the fourth oscillation, when amplitudes become smaller. As expected, single cell 
behavior differs from the average only in amplitude, significantly so for 2a deviations from 
average. But periodicity is maintained for 8 hours after IKK activation, and thus, insofar as 
the data of Nelson et al. (2004) show significant early variation in oscillation period from cell 
to cell, there have to be other sources of variability in a model with strong linear transcription. 

On the contrary, for quadratic transcription, fluctuations in the initial number of NF-kB 
molecules, affect both amplitude and periodicity of oscillations, not just their amplitude as 
in the case of linear transcription. The result is that contrary to what happens in Figure 6 
for linear transcription, where oscillations remain clear and distinct, for an initial Gaussian 
distribution of NF-kB over cells (at a fixed number of activated IKK), here oscillations in the 
number of nuclear NF-kB oscillations, are very much damped out as shown in Figure 7, as a 
result of cell-to-cell variability in both amplitude and period of oscillations. 

3.2.2 IKK fluctuations. 

Within the model parameters considered, we turn to studying fluctuations in the initial number 
of IKK which activates the cellular oscillatory response. Whatever the type of transcription, 
linear or quadratic, fluctuations in the initial amount of activated IKK lead to strong cell-to- 
cell variability. We show here our results for the case of linear transcription. 

Figure 8 shows, for fixed NF-kB = 30000, the average over many cells for a Gaussian 
distribution of the initial number of IKK (average=30000, a = 5000), and single cell response 
for IKK=25000 and IKK=35000. Several points can be made here: 

- when the initial number of IKK activated by TNFa fluctuates, not only does the ampli- 
tude of NF-kB oscillations vary, but also their periodicity 

- as a result, single cell fluctuations are washed out in the average. Average response over 
the cell population represents poorly single cell response beyond the first oscillation 

- the difference between maxima and minima of the oscillations gets larger as IKK de- 
creases from its average, while at the same time oscillatory period increases. Possible impli- 
cations of the latter result are given in the discussion. 

It is clear from the figures that whatever the source of cellular fluctuations, cell-to-cell vari- 
ability typically increases with time after the onset of IKK activation. This can be illustrated 
further by computing the standard deviation of the noise in the number of nuclear NF-kB. In 
Figure 9 we compare standard deviations as a function of time for the case of Figure 8, when 
IKK fluctuates according to a Gaussian of average and width of respectively equal 30000 and 
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5000 molecules, at fixed initial concentration of NF-kB of 30000 molecules (0.05 /j,M). On 
one hand there is af ot = < Nn 2 > — < Nn > , the total variance of nuclear NF-kB oscil- 
lations, on the other hand there is af nt = < Nn 2 > — < Nn > 2 , the intrinsic variance. Here 
Nn denotes nuclear NF-kB, the brackets denote average over intrinsic fluctuations, and the 
overline average over fluctuations in the external parameter (Elowitz et al., 2002; Swain et al., 
2002). We observe (Fig. 9) that a tot , which includes extrinsic fluctuations due to fluctuations 
in IKK, is much larger than a int except for large times (13 hours) after IKK activation. Thus 
for most of the time after IKK activation, extrinsic fluctuations dominate. The magnitude of 
G to t shortly after IKK activation reflects the significance of IKK fluctuations in determining 
cell-to-cell variability in terms of both amplitude and period of NF-kB oscillations. 

In this subsection, we have studied the effect of varying the initial amount of activated 
IKK around an average value, which we denote by IKK , while keeping the number of 
NF-kB molecules fixed. Let us now comment on how average oscillations change when 
IKK changes significantly compared to the value considered up to now. When IKK de- 
creases by 50% compared to the value chosen by Hoffmann et al. (2002), oscillations for both 
linear and quadratic transcription remain strong with an average increase in period of 20 to 30 
percent; for an increase of 50% in IKK , the period of oscillations decreases for quadratic 
transcription, whereas for linear transcription oscillations are significantly damped beyond the 
first one. In the experiments discussed (Hoffmann et al., 2002; Nelson et al., 2004) there are 
no results for the effect of varying the amount of the external stimulus, TNFa, on NF-kB 
oscillations. We note that the utility of our predictions for the average oscillations as a func- 
tion of activated IKK depends on knowing how the amount of activated IKK is related to 
the amount of external stimulant. 

4 Discussion. 

The data of Hoffmann et al. (2002) and Nelson et al. (2004) on NF-kB oscillations are 
seemingly contradictory. Hoffmann et al. (2002, suppl. material, fig. SI) observe well- 
defined oscillations for a population of cells, whereas Nelson et al. (2004) observe (for the 
identical stimulus) large cell-to-cell variability, in both amplitude and period of oscillations, 
which would normally lead to washed out average oscillations. However both groups use 
different cell lines, and experimental methods differ (for a discussion of these see the technical 
comment of Barken et al. (2005) and the response by Nelson et al. (2005)). 

From a modeling point of view, for the class of reactions considered here (cf. Table 1), we 
summarize our key results: 

- intrinsic fluctuations are small when transcription is strong (which is the case of the model of 
Hoffmann et al.(2002)), whether transcription is linear or quadratic in the number of NF-kB 
molecules. However intrinsic fluctuations can be strong, as significant as any extrinsic noise, 
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if the rate of promoter binding is small, and rate constants for transcription and translation of 
IkB become very different from the ones used by Hoffmann et al. (2002); if this is the case 
intrinsic noise can account for the observed cell-to-cell variability. Hence, it is important to 
determine transcription rate constants experimentally. 

In all cases considered, extrinsic fluctuations can be significant: 

- in the model of Hoffmann et al. (2002), fluctuations in the amount of activated IKK 
molecules can account for the observed cell-to-cell variability (Nelson et al., 2004), because 
they affect both amplitude and period of NF-kB oscillations, whatever the mode of transcrip- 
tion, linear or quadratic in the number of NF-kB molecules, of IkBql. 

- as to fluctuations in the amount of NF-kB molecules, their impact depends on the type 
of transcription, linear or quadratic. For linear transcription the amplitude only of oscillations 
varies, not their periodicity, whereas for quadratic transcription, both amplitude and periodic- 
ity vary. Only in the latter case would these fluctuations be a possible source of the observed 
cell-to-cell variability (Nelson et al., 2004) 

Fluctuations in individual cell behavior thus provide a source of damping for NF-kB 
population oscillations that is different from the one considered by Hoffmann et al. (2002), 
which depends on the existence of IkB isoforms besides IkB a. 

Our results on cell-to-cell variability of NF-kB oscillations (in the context of the model 
of Hoffmann et al., 2002) help clarify the experimental issues that need to be resolved for 
determining its origins. These requirements can be formulated as questions: what is the dis- 
tribution over the cell population of the amount of NF-kB, and of activated IKK; what are 
the rates of transcription and translation; is transcription linear or quadratic (in the right range) 
in the number of NF-kB molecules? When answers to these questions will be forthcoming, 
one will be able to determine, within a given model, the sources of cell-to-cell variability of 
NF-kB oscillations, and eventually reconcile the well-defined average oscillations observed 
by Hoffmann et al. (2002) and the large cell-to-cell variability of Nelson et al. (2004) as 
a matter of differences in the amount of fluctuations due to cell type, cell environment or 
experimental method. 

The key biological issue is the functional role of NF-kB oscillations, which is not known. 
Since NF-kB participates in regulating the expression of many genes, its oscillations can be 
transferred to the corresponding proteins. One can speculate that the number of NF-kB 
oscillations in a given time period could play the role of a threshold for setting into motion 
many aspects of cellular response. Consider Figure 8: as the level of IKK goes down the 
period of oscillations lengthens. Assume that four oscillations in the first 400 minutes after 
IKK activation represent the threshold below which signaling based on NF-kB oscillations 
becomes ineffective. In Figure 8 , the case 7^^=25000 is marginal (NF-kB is fixed). Thus 
as the amount of IKK goes lower, say to IKK =20000, and if this amount represents the 
average amount of IKK activated by an outside agent, only fluctuations in IKK will enable 
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some cells to reach the threshold of four oscillations in the required time, and thus respond 
to the external stimulation. This effect of stochasticity would lead to cellular response even 
for small, below threshold, doses of external agent, as might happen for immune response to 
a pathogen. 
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Table 1 : Reactions and Rate constants 



Rate constant 


Reaction 


Value 




7c + A^c -> 7ATc 


30/xM -1 mm -1 




iTVc 7c + Nc 


0.03mm -1 


a 7 


+ JJVc -> /ATcp 


ll.l//M -1 mm -1 


di 


JATcp -> + JJVc 


0.075mm -1 


r 4 


JA^cp -> /Kir + A^c 


1.221mm -1 


ai 


IKK + 7 c — >• 7cpr 


1.35//M -1 mm -1 


di 


7cpr — >■ IKK + 7c 


0.075mm -1 


04 


7cpr + A^c — > INcp 


30//M -1 mm -1 


C?4 


INcp — > 7cpr + A^c 


0.03mm -1 


ri 


7cpr — > IKK 


0.2442mm -1 




Ic-> 


0.00675mm -1 


tpi 


7c — * 7n 


0.018mm -1 




7n — >• 7c 


0.012mm -1 


fcl 


Nc -> iVn 


5.4mm -1 


fcbi 


iVn -»• Nc 


0.0048mm- 1 


A; 2 


INn -> JJVc 


0.8294mm- 1 


tr 2 


A^n — > A^n + 7ctr 


0.0582mm- 1 


tr 3 


7ctr -> 


0.0168mm- 1 


tri 


7ctr — > Ictr + 7c 


0.2448mm- 1 


04 


In + Nn^ INn 




d 4 


INn -> In + Nn 


0.03mm- 1 


^02 


IKK -»• 


0.0072mm- 1 



Rate constants are from Nelson et al. (2004), based on the model of Hoffmann et al. (2002). 
Notation: 

Nn= nuclear NF-kB, Nc= cytoplasmic NF-kB 
In= nuclear IkBol, Ic=cytoplasmic IkBol 
INc= cytoplasmic NF-nB-InBa complex 
INn= nuclear NF-nB-InBa complex 
INcp= complex of IKK and INc 
Icpr= complex of IKK and InBa 
Ictr= 7/t5« mRNA 
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6 Figure captions 



Figure 1. Nuclear NF-kB as a function of time for the model with linear transcription. IKK 
activation sets in after 3000 minutes. Initial total amount of NF-kB and amount of IKK are 
fixed at 0.05 fiM. The full curve is average behavior, the broken curves show three randomly 
chosen single cell oscillations. 

Figure 2. Nuclear NF-kB as a function of time for the model with linear transcription. 
Here the values of rate constants tr2 for transcription and trl for translation are 1/10 and 10 
times respectively the values given in the Table, which are used for Figure 1 (see text). IKK 
activation sets in after 3000 minutes. Initial total amount of NF-kB and amount of IKK are 
fixed at 0.05 fiM. The full and two broken curves correspond to three randomly chosen single 
cell oscillations and show the increased cell-to-cell variability as compared to Figure 1. 

Figure 3. Nuclear NF-kB as a function of time for the model variant described in section 
3.1.1. Here kb = 1200mm -1 . IKK activation sets in after 3000 minutes. Initial total amount 
of NF-kB and amount of IKK are fixed at 0.05 /iM. The full curve is the average behavior 
of NF-kB, the broken curve represents the standard deviation due to intrinsic noise. 

Figure 4. Nuclear NF-kB as a function of time for the model variant described in section 
3.1.1. Here kb = \2min~ 1 ; transcription is weak. The full curve is average cell behavior, the 
broken line is the standard deviation a. 

Figure 5. Nuclear NF-kB as a function of time for the model variant described in section 
3.1.1. Here kb = 12mm -1 ; transcription is weak. The curves correspond to the behavior of 
three single cells. Oscillations differ in both amplitude and phase. 

Figure 6. Nuclear NF-kB as a function of time for the model with linear transcription. 
IKK activation sets in after 3000 minutes. The amount of initially activated IKK is 0.05 fiM. 
The full curve represents an average over 1000 single cells for an initial amount of NF-kB 
following a gaussian distribution of average number equal to 30000 (corresponding to 0.05 
fiM for the cell volume chosen), and width a = 5000. The broken curves correspond to 
one single cell behavior for, by decreasing amplitude, NF-kB= 35000, 25000, and 20000 
molecules. 

Figure 7. Nuclear NF-kB as a function of time for the model with quadratic transcription. 
IKK activation sets in after 3000 minutes. The amount of initial activated IKK is 0.05 fj,M. 
The full curve represents an average over 1000 single cells for an initial amount of NF-kB 
following a Gaussian distribution of average number equal to 30000 (corresponding to 0.05 
fiM for the cell volume chosen), and width a = 5000. 

Figure 8. Nuclear NF-kB as a function of time for the model with linear transcription. 
IKK activation sets in after 3000 minutes. The amount of initial NF-kB is fixed at 0.05 [iM. 
The full curve represents an average over 1000 single cells for an initial amount of IKK 
following a Gaussian distribution of average number equal to 30000 (corresponding to 0.05 
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/jlM for the cell volume chosen), and width a = 5000. The broken curves correspond to single 
cell behavior for, by decreasing amplitude, IKK= 25000 and 35000 molecules. 

Figure 9. Standard deviations (in /iM) of nuclear NF-kB oscillations for the case consid- 
ered in Figure 3, as a function of time. IKK activation sets in after 3000 minutes. a tot (broken 
line) represents the total noise, a int (full line) the intrinsic noise (see text). 
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